OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Exome Sequencing Analysis Identifies Compound 
Heterozygous Mutation in ABCA4 in a Chinese Family 
with Stargardt Disease 

Yu Zhou''^'*', Siyu Tao^', Hui Chen" '^ Lulin Huang^'^'^ Xiong Zhu^'^'^ Youping Li^ Zhili Wang^ 
He Lin^ '^ Fang Hao^ Zhenglin Yang^'^'^'^ Liya Wang^*, Xianjun Zhu^'^'^'^'^* 

1 Sichuan Provincial Key Laboratory for Human Disease Gene Study and Institute of Laboratory Medicine, Sichuan Academy of Medical Sciences and Sichuan Provincial 
People's Hospital, Chengdu, Sichuan, China, 2 School of Medicine, University of Electronic Science and Technology of China, Chengdu, Sichuan, China, 3 Henan Eye 
Hospital and Henan Eye Institute, People's Hospital of Zhengzhou University, Zhengzhou, Henan, China, 4 Department of Ophthalmology, Sichuan Academy of Medical 
Sciences and Sichuan Provincial People's Hospital, Chengdu, Sichuan, China, 5 Laboratory Animal Institute, Sichuan Academy of Medical Sciences and Sichuan Provincial 
People's Hospital, Chengdu, Sichuan, China, 6 Chengdu Institute of Biology, Chinese Academy of Sciences and Chinese Academy of Sciences Sichuan Translational 
Medicine Research Hospital, Chengdu, Sichuan, China, 7 Key Laboratory for Neurolnformation of Ministry of Education, University of Electronic Science and Technology of 
China, Chengdu, Sichuan, China, 8 College of Life Sciences and Engineering, Xinan Jiaotong University, Chengdu, Sichuan, China 



Abstract 

Stargardt disease is the most common cause of juvenile macular dystrophy. Five subjects from a two-generation Chinese 
family with Stargardt disease are reported in this study. All family members underwent complete ophthalmologic 
examinations. Patients of the family initiated the disease during childhood, developing progressively impaired central vision 
and bilateral atrophic macular lesions In the retinal pigmental epithelium (RPE) that resembled a "beaten-bronze" 
appearance. Peripheral venous blood was obtained from all patients and their family members for genetic analysis. Exome 
sequencing was used to analyze the exome of two patients 111, 112. A total of 50709 variations shared by the two patients 
were subjected to several filtering steps against existing variation databases. Identified variations were verified In all family 
members by PCR and Sanger sequencing. Compound heterozygous variants p.Y808X and p.G607R of the ATP-blndIng 
cassette, sub-family A (ABCl), member 4 (ABCA4) gene, which encodes the ABCA4 protein, a member of the ATP-blndIng 
cassette (ABC) transport superfamlly, were Identified as causative mutations for Stargardt disease of this family. Our findings 
provide one novel ABCA4 mutation in Chinese patients with Stargardt disease. 
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Introduction 

Stargardt disease (STGD), which is also known as juvenile 
macular degeneration, was first reported by Karl Stargardt in 
1909. It is one of the most common hereditary retinal dystrophies 
with an estimated prevalence of at least 1:10,000 [1,2,3]. It 
presents with a progressive and significant loss in central vision in 
the first or second decade of life. However, fundus examination is 
frequently normal early in the course of disease, even when 
patients already complain of vision loss. At this stage, the clinical 
diagnosis of Stargardt disease may be missed. Later on, typical 
fundus manifestations arise, including pigment mottiing, frank 
macular atrophy, a buU's eye maculopathy and fundus flecks in the 
macular and the perimacular region [4]. However, it should be 
noted that Stargardt disease presents with highly variable 
phenotypes. Histologically, Stargardt disease is associated with 
significant loss in photoreceptor cells and a massive deposition of 



lipofuscin-like material in the retinal pigment epithelium, which 
has also been observed in aging human eyes [6,7]. 

STGD is predominantiy inherited as an autosomal recessive 
trait, although an autosomal dominant form has been also 
described [5]. Both sexes are equally affected. The STGD gene 
has been mapped to the short arm of chromosome 1 [8] in a 
narrow genetic interval, subsequentiy assigned to band p22.1 [9], 
now known as ATP-binding cassette, sub-family A (ABCl), 
member 4 (ABCA4) [10,1 1]. The gene for an autosomal dominant 
disorder with a similar phenotype has been reported on 
chromosome 6 [12]. Autosomal dominant families linked to this 
locus are at least 50 times less common than families that are 
consistent with autosomal recessive inheritance [13,14]. Reces- 
sively inherited Stargardt disease is likely to be monogenic. Rare 
cases of STGD or "Stargardt-like" disease phenotypes have been 
reported with mutations in PRPH2 [15,16], VMD2 [17], EL0VL4 
[5,18,19] and PROMl [20]. These genes, as well as ABCA4, are 
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Figure 1. Pedigree of family 2048 witli Stargardt. Solid symbols 
indicated affected individuals. Open symbols indicated unaffected 
individuals and arrow indicates the proband. 
doi:1 0.1 371 /journal.pone.0091 962.g001 

also associated with clinically distinct phenotypes including 
retinitis pigmentosa, cone/rod dystrophy and pattern dystrophy. 

Known candidate genes for Stargardt disease such as "ABCA4" 
contain many exons and there are hundreds of identified 
mutations. The cost and time requirement for mutation screening 
of all coding exons by Sanger sequencing would equal or exceed 
that of high-throughput next generation sequencing (NGS) 
analysis. In this study, we applied next generation sequencing 
technology to identify the disease-causing gene in this family as 
part of a large cohort study for retinal diseases. Next-generation 
sequencing, in particular whole-exome sequencing (WES) can now 
be performed rapidly and at minimal cost, allowing analysis of the 
coding regions (exome) of the human genome in single individuals 
or small families, including patients in whom a clear genotype- 
phenotype correlation is absent or for clinically and genetically 
heterogeneous conditions [21,22,23]. 

In the present study, disease-associated mutations were identi- 
fied by WES of the two affected siblings followed by vahdation in 
the family affected by Stargardt disease. Our results identified two 
compound heterozygous disease-segregating mutations, 
C.C2424G, P.Y808X and C.G1819A, p.G607R, in the ABCA4 
gene. To exclude the possibihty that these mutations were 
polymorphisms, DNA samples of 1000 unaffected individuals 
were also screened for these mutations. 

Materials and Methods 

Subjects and Clinical Assessment 

Study approval was obtained from the Institutional Review 
Boards of Sichuan Academy of Medical Sciences and Sichuan 
Provincial People's Hospital and Henan Eye Hospital and Henan 
Eye Institute, People's Hospital of Zhengzhou University. Written 
informed consent was obtained in accordance with the 



Declaration of Helsinki for all subjects enrolled in this study. For 
minors, written consent was obtained from the father. Five family 
members were evaluated by a retina specialist and underwent 
complete ophthalmological assessment that included visual acuity 
measurement, fundus photography, fundus fluorescein angiogra- 
phy (FFA), multifocal electroretinogram (mfERG), optical-coher- 
ence tomography (OCT) and computerized visual field testing. In 
the 1000 normal matched controls, aU individuals underwent an 
eye examination and no signs of eye diseases were observed. 
Venous blood samples were obtained from all subjects in EDTA 
Vacutainers. 

DNA Extraction 

All genomic DNA was extracted from peripheral blood using a 
blood DNA extraction kit according to the protocol provided by 
the manufacturer (TianGen, Beijing, China). DNA samples were 
stored at — 20°C until used. DNA integrity was evaluated by 1% 
agarose gel electrophoresis. 

Exome Sequencing 

Exome sequencing was employed in this study to identify the 
disease-associated genes. Exome sequencing was performed on 
DNA samples of the two patients (family members II 1 and 112) by 
Axeq Technology Inc., Seoul, Korea. Each sequenced sample was 
prepared according to the lUumina protocols. Briefly, one 
microgram of genomic DNA was fragmented by nebulization, 
the fragmented DNA was repaired, an 'A' was ligated to the 3' 
end, lUumina adapters were then ligated to the fragments, and the 
sample was size selected aiming for a 350 — 400 base pair product. 
The size-selected product was PGR amplified, and the final 
product was validated using the Agilent Bioanalyzer. Streptavidin 
beads were used to capture probes containing the targeted regions 
of interest; non-specific binding was then washed out. Then, the 
sequencing libraries were enriched for the desired target using the 
lUumina Exome Enrichment protocol and the enriched library 
validation for quality control analysis were performed by the Axeq 
Technology. For clustering and sequencing, genomic DNA 
lUumina TruSeq Exome Capture System (62 Mb) was used to 
coUect the protein coding regions of human genome DNA. It 
covered 20794 genes and 2011 21 exons in the Consensus Coding 
Sequence Region database, approximately 97.2% of CCDS exons 
or 96.4% of RefSeq exons were captured. (http://www.illumina. 
com/apphcations/sequencing/ targeted resequencing.Umn) Each 
captured library was then loaded onto the lUumina Hiseq2000 
sequencer, and we performed high-throughput sequencing for 
each captured library to ensure that each sample met the desired 
average sequencing depth. Raw image files were processed by 



Table 1. Phenotype and genotype of the family members. 





Family 
Member 


Age 


Gender 


Disease 

Duration 

(years) 


Visual Acuity 
(OD, OS) 


IVIutatlon(s) 




Original reports 
described 












Nucleotide change 


Effect 




11 


43 


M 


0 


20/20, 20/20 


18190A 


Gly 607 Arg 


Andrea Rivera 


12 


43 


F 


0 


20/20, 20/20 


24240G 


Tyr 808* 


N/A 


111 


13 


M 


3 


20/400, 20/400 


18190A/24240G 


Gly 607 Arg,/Tyr 808* 


Andrea Rivera, N/A 


112 


18 


F 


8 


20/400, 20/400 


18190A/24240G 


Gly 607 Arg,/Tyr 808* 


Andrea Rivera,N/A 


113 


20 


F 


0 


20/20, 20/20 


1819G>A 


Gly 607 Arg 


Andrea Rivera 



N/A, Not available; M, Male; F, Female. 
doi:l 0.1 371 /journal.pone.0091 962.t001 
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(Relative To Internal Normal) 
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nV/dog' 




Ring Amp.P1 Amp.P1 Amp.N1 PeT.P1 
(nV/defl'l [pV] [pV] [msl 



56 0.48 0.04 47.1 

18.6 0.22 0.01 50.0 

21.5 0.34 0.05 46.1 

0.55 0.01 0.01 19.6 

3.67 0.10 0.01 49.0 

3.12 0.11 0.12 49.0 



3D Amplitudes P1 (b) 
(Relative To Internal Normal) 



Amplitudes P1(b) 



5D0f V7?3f 




(ISO 




20nis/di\' 



Ring Amp.P1 Amp.P1 Amp.N1 PeT.P1 
(nV/degn [pV] [pV] rmsl 



50.8 
23.6 
8.04 
4.42 
5.19 
3.31 



0.437 0.175 37.2 
0.273 0.119 58.8 
0.129 0.106 50.0 
0.095 0.062 38.2 
0.144 0.036 45.1 
0.121 0 088 40 2 
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Figure 2. Representative photograpKis of patients of family 2048. (A) Fundus photographs showing pigment mottling and yellow-white 
flecks in both maculae. (B) Fluorescein angiography (FA) images showing the hyperfluorescent flecks extended to the midperipheral retina and 
fluorescence blocking by the pigment mottling in the mcular. (C) mfERG records showing severe depressed central waveform and significant 
paracentral/pereferral loss of retinal response. (D) Macular OCTs showing hyper-reflective deposits within the RPE layer and the level of the outer 
segments of the photoreceptors, thinning of the retinal outer layers and enhanced choroidal reflectivity associated with overlying atrophic retina. 
doi:1 0.1 371 /journal.pone.0091 962.g002 



Illumiiia base calling Software 1.7 for base calling with default 
parameters and the sequences of each individual were generated 
as 90-bp pair-end reads. 

Reads, Mapping, and Variant Detection 

The high-quality sequencing reads were aligned to the human 
reference genome (NCBI build 37.1/ hgl9) with SOAPaligner 
(soap2.21). Based on the SOAP alignment results, SOAPsnp vl.05 
was used to assemble the consensus sequence and call genotypes in 
target regions. Data were provided as hsts of sequence variants 
(SNPs and short Indels). For SNP quality control, we filtered 
SOAPsnp results as follows: (i) Base quality is more than 20; (ii) 
Depth is between 4 and 200; (iii) Estimate copy number is equal or 
less than 2; (iv) The distance between two SNPs must be longer 
than 4. Small Indel detection was performed using the 
UnifiedGenotyper tool from GATK (version vl. 0.4705) after all 
the high-quality reads were aligned to the human reference 
genome using BWA (version 0.5.9-rl6). SNP and Indel detection 
were performed only on the targeted exome regions and flanking 
regions within 200 bp. 

Filtering and Annotation 

The detected variants were annotated and filtered based on four 
databases; i.e., NCBI CCDS (http:/ /www.ncbi.nlm.nih.gov/ 
CCDS/CcdsBrowse.cgi), RefSeq (http:/ /www.ncbi.nlm.nih.- 
gov/RefSeq/), Ensembl (http://www.ensembl.org), and Encode 
(http://genome.ucsc.edu/ENCODE). Four major steps were 
taken to prioritize all the high-quality variants: (i) variants within 
intergenic, intronic, and UTR regions and synonymous mutations 
were excluded from downstream analysis; (ii) variants in 
dbSNP137(http://www.ncbi.nlm.nih.gov/projects/SNP/), 1000 
Genome project (ftp://ftp.1000genomes.ebi.ac.uk/voll/ftp), YH 
Database (http://yh.genomics.org.cn/), HapMap Project (ftp:// 
ftp.ncbi.nlm.nih.gov/hapmap) and our in-house database, which 
was generated by our laboratory using 1600 whole exome 
sequencing data, were excluded; (iii) possible damaging impacts 
of each variant on protein structure /function were predicted by 
SIFT (http://sift.bii.astar.edu.sg/) and Polyphen2 (http:// 
genetics.bwh.harvard.edu/pph2/); (iv) gene Ontology (http:// 
www.geneontology.org) and KEGG Pathway Annotations 



(http://www.ebi.ac.uk./clustalw) were used to predict the biolog- 
ical function of each putative gene. 

Variants Validation 

After filtering against multiple databases, Sanger sequencing 
was used to determine whether any of the remaining variants co- 
segregated with the disease phenotype in this family. Primers 
flanking the candidate loci were designed based on genomic 
sequences of Human Genome (hgl9/build37.1) and synthesized 
by Invitrogen, Shanghai, China: ABCA4-exonl3-F, tgagttccgagt- 
caccctgt; ABCA4-exonl3-R, gtcagagctccatgctctcc; ABCA4-ex- 
onl6-F, ctctacctcgagggcatctg; ABCA4-exonl6-R, ggctggggatct- 
gaagaact. Genotyping for C.C2424G and C.G1819A in the 
family members was then confirmed by direct polymerase chain 
reaction (PGR) and analyzed on an ABI 3730XL Genetic 
Analyzer. Sequencing data were compared with the Human 
Genome database. 

Results 

Clinical Presentation of Family 2084 

A two-generation family (family 2084) from Henan Province of 
China was recruited in this study (Figure 1). Ophthalmic 
examinations identified two affected individuals as Stargardt 
disease patients among the five examined family members. 
Affected members of this family exhibited similar clinical features. 
The two affected siblings presented with an early-onset markedly 
decreased vision acuity (OD: 20/400, OS: 20/400) in both eyes 
and an increasing difficulty to adapt in the dark (Table 1). Fundus 
examination showed some pigment mottling and yellow-white 
flecks in both maculae, normal caliber of the retinal vessels, but no 
pigmented bone spicules in the retinal periphery (Figure 2A). The 
fluorescein angiogram displayed hyperfluorescent flecks, which 
extended to the midperipheral retina and fluorescence blocks 
formed by the pigment mottling in the macula (Figure 2B). 
Multifocal Electroretinogram showed severe depressed central 
waveform and significant loss of paracentral/peripheral retinal 
response (Figure 2C). The macular OCT showed hyper-reflective 
deposits within the RPE layer and the level of the outer segments 
of the photoreceptors, thinning of the retinal outer layers and 
enhanced choroidal reflectivity (Figure 2D). 



Table 2. Number of candidate SNP/lndels filtered against several public variation databases and the in-house data. 







Feature_SNP shared by Case 111 
and Case 112 


Featurejndel shared by Case 111 
and Case 112 


TotaLnumber 


49535 


1174 


FunctlonaLSNP/lndels 


5843 


334 


Filtered_DBsnp1 35common/indel 


754 


228 


Filtered_DBsnp/indel_1 000gene(201 1 ) 


753 


226 


Filtered_DBsnp_1 OOOgene_Hapmap 


752 


226 


F[ltered_DBsnp_1000gene_Hapmap_YH 


748 


226 


Filtered in-House Data 


25 


0 



doi:1 0.1 371 /journal.pone.0091 962.t002 
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Table 3. Candidate exome sequence variants shared by all affected individuals and filtered by database. 



Chr:Position Gene Mutation 



Mode 



Maternal Paternal SIFT/ 
Mutation type dbSNP ID Allele Allele Polyphen2 



Stargardt 
Gene 



chr01:94520830 ABCA4 exonl6:c.C2424G:p.Y808X heterozygous stopgain Novel Yes 

chr01:9452825 ABCA4 exon13:c.G1819A:p.G607R heterozygous nonsynonymous rs61749412 



Yes 



Stop codon, 
N/A 

Damaging/ 

probably 

damaging 



N/A 



N/A, Not available. 

doi:l 0.1 371 /journal.pone.0091 962.t003 



Whole-exome Sequencing 

By exome sequencing of patient 11:1 and patient II;2, we 
generated 6.9 and 7.2 billion bases of sequence with average 
throughput depth of target regions Ul.Gx and 116. 4x, respec- 
tively. Approximately 99.8% and 99.8% of initial mappable reads 
were able to pass our thresholds for calling SNPs and short 
insertions or deletions (Indels). The mean read depth of target 
regions was 52.3 x and 54.7 x, respectively. In 11:1 and 11:2, we 
separately identified 20563 and 20641 SNPs in the coding regions 
(9377 and 9407 nonsynonymous SNPs, 10663 and 10689 
synonymous SNPs, 523 and 545 other types of SNPs, respectively), 
423 and 413 coding Indels that may affect amino acid sequence, 
respectively. 

After identification of variants, we identified 5843 functional 
SNP and 334 functional Indels that were shared by these two 
patients (Table 2). We focused only on the functional SNP/Indel, 
including non-synonymous variants (NS), splice acceptor and 
donor site mutations (SS), and frameshift coding-region insertions 
or deletions (Indels), which were more likely to be pathogenic than 
others, especially those in homozygous or multiple heterozygous 
status. We then compared these variants in two affected members 
with the dbSNP135, 1000 Genome Project, HapMap project, YH 
database and our in-house database, which generated by our 



Patient 11:1 



Patient 11:2 



Father 1:1 



Mother 1:2 



\AA./waa/1AaA hkfMVUh 

GAGTACCTGGTT GGGGCGGGTTTG 

\AAaAA.aaAAaA .MaLa/Iaaa/I 



G AGTACCTG GT T 



GGGGCGGGTTTG 



GAGTACCTGGTT GGGGCGGGTTTG 




t m 




GAGTACCTGGTT GGGGCGGGTTTG 



2424 OG 



1819 G>A 



Figure 3. Mutation identification of ABCA4 gene. Electrophero- 
gram analysis of ABCA4 in family 2048 showing the compound 
heterozygous mutations (C.C2424G and C.G1819A) co-segregated with 
the phenotype. Ill and 112 patients harbored compound heterozygous 
C.C2424G and C.G1819A mutations of the ABCA4 gene. C.C2424G 
mutation was carried by the mother 12 while C.G1819A mutation was 
carried by the father 11. 
doi:1 0.1 371/journal.pone.0091 962.g003 



laboratory using 1600 whole exome sequencing data (table 2). The 
in-house data include whole exome variants data from people 
without any eye disease, therefore we can exclude the variants with 
high frequency in people without any eye disease. 

Under the autosomal recessive model, variants satisfying a 
recessive homozygous inheritance model were not identified. This 
led us to investigate the possibility of recessive compound 
heterozygous inheritance. Using this model, the filtered data was 
narrowed to 25 heterozygous variants. We then compared these 
variants with reported retina genes (https:/ /sph.uth.edu/Retnet/). 
In both patients we found two mutations C.C2424G (p.Y808X) 
and C.G1819A (p.G607R) satisfying a recessive compound 
heterozygous inheritance model (Table 3) in the gene ABCA4 
(NM_000350.2). When we checked the human gene mutation 
database (http:/ /www.hgmd.org/), we found that mutation 
p.G607R was reported by Andrea Rivera in 2000 [24]. Sanger 
sequencing confirmed these two mutations in the two affected 
siblings and demonstrated that their parents were unaffected 
carriers of Y808X (father) and G607R (mother) mutations, 
showing complete co-segregation of the mutations with the disease 
phenotype (Figure 3). The two mutations described above were 
absent in 1000 ethnicity-matched control samples screened by 
direct sequencing. These data, together with the clinical presen- 
tation of the two affected siblings, demonstrated that p.G607R and 
p.Y808X variants in the gene ABCA4 was responsible for Stargardt 
disease in this family. 

Mutation Detection and Analysis 

SIFT was used to predict how the identified amino acid 
substitutions would affect protein function. The previously 
reported mutation p.G607R of ABCA4 was predicted to be 
damaging and the novel mutation C.C2424G, p.Y808X in the 
affected families introduced a stop codon, which removed 1465 
amino acids from the ABCA4 protein (2273 amino acids), 
according to GenBank accession number NM_000350.2. There- 
fore, this novel mutation is likely a nuU allele. 

Polypheii2 was used to explore sequence homology and the 
physical properties of corresponding affected amino acids. 
Mutation p.G607R, located within exon 13, results in changes 
of a hydrophHic glycine to an arginine at position 607, which may 
lead to a damaging replacement with a score of 0.99 (sensitivi- 
ty:0.72, specificity:0.97) Figure 4A. Mutation Y808X, located 
within exon 16, results in a nonsense mutation, and the mRNA 
with a premature stop codon is likely to be degenerated by the 
nonsense-mediated mRNA decay response, thus leading to a 
decrease in ABCA4 expression. As shown in Figure 4A, both amino 
acid changes affect highly conserved residues. 

We then used TMHMM2.0 to predict the ABCA4 protein 
structure. Our result showed that the protein was organized in two 
structurally related tandem-arranged halves with each half 
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Mutated 

H. sapiens 

D . rerio 

G.gallus 

M.mus cuius 

R . norveglcus 

M.mulatta 

P . troglodytes 

C . lupus 

B . taurus 



GFAYLQDMVEQRIT- 
GFAYLQDMVEQGIT- 
GFAYIiQDLIEHGIL- 
GFAYLQDMIEHGI I- 
GFAYLQDMVEQGIV- 
GFAYIiQDMVEQGIV- 
GFAYLQDMVEQGIT- 
GFAYLQDMVEQGIT- 
GFAYLQDMIEQGI T- 
GFAYLQDMVEHGIT- 



RSQVQAEAPVGI 
RSQVQAEAPVGI 
KLHTGHDWPLGV 
KTQTNTDVALGI 
KSQMQAEPPIGV 
RSQTPAEPPIGV 
RSQVQAEAPVGI 
RSQVQAEAPVGI 
RSQAQVEVPVGI 
RSQAQEEVPVGI 



LSPVAFGFGTEXXXXXXXXXXXXXXXXX 
LSPVAFGFGTEYLVRFEEQGLGLQWSNI 
LSPVAFGFGTEYLSRYEEQGLGLQWDNI 
LSQVAFGFGTEYLSRYEEQGLGLQWGNI 
LSSVAFGFGTEYLVRFEEQGLGLQWSNI 
LSPVAFGFGTEYLVRFEEQGLGLQWSNI 
LSPVAFGFGTEYLVRFEEQGLGLQWSNI 
LSPVAFGFGTEYLVRFEEQGLGLQWSNI 
LSPVAFGFGTEYLARFEEQGLGLQWSNI 
LSPVAFGFGTEYLAXFEEQGVGLQWSNI 



G607R 



Y808X 



B 



Disc 
Lumen 




Cytoplasm 



Figure 4. Topological organization of ABCA4 and conservation analysis. (A) Protein alignment showed conservation of residues ABCA4 
Y808X and G607R across nine species. These two mutations occured at an evolutionarily conserved amino acid. The ABCA4 mutation G607R occured 
in the NBD1 while Y808X occured in the disc lumen region between TIV1D1. (B) Topological organization of ABCA4 in the disc membrane was shown. 
The domain organization included exocytoplamic domain (ECD), nucleotide binding domain (NBD) and transmembrane domain (TMD). 
doi:10.1371/journal.pone.0091962.g004 



containing transmembrane domains (TMD) followed by nucleo- 
tide binding domains (NBD). Our result also indicated that G607R 
occurs in the ECDl domain, while Y808X occurs in the disc 
lumen region between TMDl (Figure 4B). 

Discussion 

This study identifes novel compound heterozygous mutations in 
the ABCA4 gene as a cause of Stargardt's disease. STGD accounts 
for approximately 7% of all retinal dystrophies. It is one of the 
most common genetic forms of juvenile or early adult onset 
macular degeneration. This condition affects the central retina 
(macula) with a variable phenotype and a variable age of onset and 
severity. 

The ABCA4 gene, located at the chromosome lp22.1 with 50 
exons, is a large glycoprotein with 2,273 amino-acid and 
organized into two structurally related tandem-arranged halves 
with each half containing a transmembrane domain (TMD) 
followed by a nucleotide binding domain (NBD) [25,26,27]. Both 
the N and C halves are predicted to have a single membrane- 
spanning segment followed by a large exocytoplasmic (extracellu- 
lar/lumen) domain (ECD), five membrane-spanning segments and 
a nucleotide-binding domain (NBD) [27,28]. A highly conserved 



VFVNFA motif near the C-terminus has been shown to play an 
essential role in the folding of ABCA4 into a functional protein 
[29]. 

The ABCA4 protein is localized in cone and rod photoreceptor 
outer segments [30]. The normal function of ABCA4 is to facilitate 
the transport of aU-trans-retinal from the outer segment disk to the 
outer segment cytoplasm in the form of a mono-substituted 
phospholipid known as N-retinylidene- phosphatidylethanolamine 
(N-ret-PE) [31,32]. When ABCA4 is defective, N-ret-PE irrevers- 
ibly forms a toxic, insoluble, bisretinoid known as A2PE which is 
deposited in retinal pigment epithelium (RPE) cells during the 
process of disc shedding and phagocytosis, eventually leading to 
cell death and macular degeneration [33]. 

Compared to wild-type mice, Abca4 knockout mice show 
significant light-dependent changes in lipids and an accumulation 
of lipofuscin deposits in the RPE cells. Biochemical analysis of the 
lipofuscin deposits from Abca4 knockout mice show elevated levels 
of several fluorescent diretinoid compounds, including A2E, a 
diretinal pyridinium compound known to be a major component 
of lipofuscin, all-trans retinal dimer, and related diretinal 
compounds [34,35,36,37,38,39]. In addition, Abca4 knockout 
mice, like individuals with Stargardt disease, show a delay in dark 
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adaptation consistent with the delayed removal of aU-trans retinal 
from outer segments following photobleaching [31]. 

Since the first report of mutations in the ABCA4 gene by 
Allikmets in 1997, over 800 disease-causing mutations have been 
identified to date in ^5 G4 ^-associated phenotypes and more than 
half of these have been detected only once [40]. The ABCA4 
mutation spectrum includes missense, nonsense, splice-site, frame- 
shift, small deletion and insertion mutations, although the 
approximately 80% of reported changes are missense mutations. 
It is now generally believed that mutations in ABCA4 result in a 
spectrum of related retinal dystrophies, including STGD, bull's eye 
maculopathy [41,42], retinitis pigmentosa [43,44,45,46,47], cone 
rod dystrophy and age-related macular degeneration [48]. 
Numerous genetic studies on STGD patients have revealed that 
the disease-associated ABCA4 alleles are extraordinarily heteroge- 
neous. It has been estimated that only half of STGD cases have 
two known or putative disease-[:ausing ABCA4 mutations on 
separate alleles, nearly one third of cases have a single mutation 
and the remaining have no definite or probable disease-causing 
ABCA4 mutations [49]. 

In our study compound heterozygous mutations p.Y808X and 
p.G607R of the ABCA4 gene were identified. The previously 
reported ABCA4 p.G607R mutation [24,50] in exon 13 is a single 
nucleotide polymorphism, rs6 17494 12, which is predicted prob- 
ably to be damaging to protein function (PolyPhen2 scores close to 

References 

1. Wclcbcr RG (1994) Stargardt's macular dystrophy. Arch Ophthalmol 112: 

752-754. 

2. Walia S, I'ishman (iA (2009) Natural history ot'phcnotypic changes in Stargardt 
macular dystrophy. Ophthalmic Genet 30: 63—68. 

3. Stargardt K (1909) Uber familiare, progressive degeenration under makulage- 
gend des augen. Albrecht von Graefes Arch Ophthalmol: 71:534—550. 

4. Franceschetti A, Francois J (1965) [Fundus flavimaculatus] . Arch Ophtalmol 
Rev Gen Ophtahnol 25: 505-530. 

5. Zhang K, Kniazcva M, Han M, Li W, Yu Z, et al. (2001) A 5-bp deletion in 
ELOVL4 is associated with two related forms of autosomal dominant macular 
dystrophy. Nat Genet 27: 89-93. 

6. Bimbach GD, Jar\'daincn M, Possin DE, Milam AH (1994) Histopatholog^' and 
immunocytochcmistry of the neurosensory retina in fundus llavimaculatus. 
Ophdiahnology 101: 1211-1219. 

7. Steinmetz RL, Gamer A, Maguire JI, Bird AG (1991) Histopathology of 
incipient fundus flavimaculatus. Ophthalmology 98: 953—956. 

8. Kaplan J, Gerber S, Larget-Ket D, RozetJM, Dollfus H, et al. (1993) A gene for 
Stargardt's disease (fundus flavimaculatus) maps to the short arm of chromosome 
1. Nat Genet 5: 308-31 1. 

9. RozetJM, Gerber S, Perrault 1, Camuzat A, Galvas P, et al. (1996) Structure 
and physical mapping of DRl, a r.\T.\-binding protein-associated phospho- 
protein gene, to chromosome lp22.1 and its exclusion in Stargardt disease 
(STGD). Genomics 36: 554-556. 

10. Nasonkin I, Illing M, Kochler MR, Schmid M, Molday RS, et al. (1998) 
Mapping of the rod photoreceptor ABC transporter (ABCR) to Ip21-p22.1 and 
identification of novel mutations in Stargardt's disease. Hum Genet 102: 21—26. 

11. Allilcmets R, Singh N, Sun H, Shroyer NF, Hutchinson A, et al. (1997) A 
photoreceptor cell-specific AfP-binding transporter gene (ABCR) is mutated in 
recessive Stargardt macular dystrophy. Nat (ienct 15: 236 — 246. 

12. Houlston RS, Chcadlc J, Dobbins SE, Tcnesa A, Jones AM, et al. (2010) Meta- 
analysis of three genome-wide association studies identifies susceptibility loci for 
colorectal cancer at lq41, 3q26.2, 12ql3.13 and 20ql3.33. Nat Genet 42: 
973-977. 

13. Stone EM, Nichols BE, Kimura AE, Weingeist TA, Drack A, et al. (1994) 
Glinical features of a Stargardt-like dominant progressive macular dystrophy 
with genetic linkage to chromosome 6q. Arch Ophthalmol 112: 765 — 772. 

14. Edwards AO, Micdziak ,\, Vrabec T, Vcrhoeven J, Acott TS, et al. (1999) 
Autosomal dominant Stargardt-like macular dystrophy: 1. Clinical character- 
ization, longitudinal follow-up, and evidence for a common ancestry in families 
linked to chromosome 6ql4. AmJ Ophthalmol 127: 426—435. 

15. Goco RM, Telleria JJ, Sanabria MR, Rodriguez-Rua E, Garcia MT (2010) 
PRPH2 (Peripherin/RDS) mutations associated with different macular dystro- 
phies in a Spanish population: a new mutation. Eur J Ophthalmol 20: 724—732. 

16. Poloschek CM, Bach M, Lagreze WA, Glaus E, Lemke JR., et al. (2010) ABCA4 
and ROMl: implications for modification of the PRPH2-associated macular 
dystrophy phenotype. Invest Ophthalmol Vis Sci 51: 4253 — 4265. 

17. Petrukhin K, Koisti MJ, BakaU B, Li 'W, Xie G, et al. (1998) Identification of die 
gene responsible for Best macular dystrophy. Nat Genet 19: 241 — 247. 



1.0). Through the analysis of membrane topology by 
TMHMM2.0, we found that the p.G607R mutation was in the 
ABCA4 ECDl region, which is involved in stacking interactions 

with the adenine ring of ATP [51]. The no\ el stopgain p.Y808X 
mutation in exon 16 was detected in a heterozygous state, close to 
the previously reported mutation p.G818E [11]. It was predicted 
to result in a truncated protein, which severely impaired ABCA4 
protein function. This novel compound heterozygous mutation 
was absent from public databases such as 1000 genomes or Exome 
Variant Server, excluding them as common polymorphisms. 

In summary, we have reported the clinical and genetic 
characteristics of a Chinese family with Stargardt disease by 
WES. To identify pathogenic variants, we analyzed these variants 
by subjecting them to an analytical pipeline for high-confidence 
variant calling, annotation and filtration and finally identified 
novel compound heterozygous mutations in ABCA4. Our study 
provides another compound heterozygous mutation in ABCA4 for 
Stargardt disease. 

Author Contributions 

Conceived and designed the experiments: YZ LW XZ ZY. Performed the 
experiments: YZ ST HC XZ YL HL FH ZLW. Analyzed the data: YZ ST 
HC HL XZ. Contributed reagents/materials/analysis tools: ST LW HC. 
Wrote the paper: YZ XZ. 



18. Edwards AO, Donoso LA, Rittcr R 3rd (2001) A novel gene for autosomal 
dominant Stargardt-like macular dystrophy with homology to the SUR4 protein 
family. Invest Ophthalmol Vis Sci 42: 2652-2663. 

19. Vasireddy V, Wong P, A^'yagari R (2010) Genetics and molecular pathology of 
Stargardt-like macular degeneration. Prog Retin Eye Res 29: 191—207. 

20. Yang Z, Chen Y, Lillo C, Ghien J, Yu Z, et al. (2008) Mutant prominin 1 found 
in patients with macular degeneration disrupts photoreceptor disk morphogen- 
esis in mice. J Clin Invest 118: 2908-2916. 

21. Nielsen R, Paul JS, Albrechtscn A, Song YS (2011) Genotype and SNP caUing 
from next-generation sec^uencing data. Nat Rev Genet 12: 443—451. 

22. Ku CS, Cooper DN, Polychronakos CI Naidoo N, Wu M, ct al. (2012) Exome 
sec^ueneing: dual role as a discover)' and diagnostic tool. .Ann Neurol 71: 5—14. 

23. Casals F, Idaghdour Y, HussinJ, Awadalla P (2012) Next-generation sequencing 
approaches for genetic mapping of complex diseases. J Neuroimmunol 248: 
10-22. 

24. Rivera A, White K, Stohr H, Steiner K, Hemmrich N, et al. (2000) A 
comprehensive survey of sequence variation in the ABCA4 (ABCR) gene in 
Stargardt disease and age-related macular degeneration. AmJ Hum Genet 67: 
800-813. 

25. Molday RS (2007) ATP-binding cassette transporter ABCA4: molecular 
properties and role in vision and macular degeneration. J Bioenerg Biomembr 
39: 507-517. 

26. Dhng M, Molday LL, Molday RS (1997) The 220-kDa rim protein of retinal rod 
outer segments is a member of the ABC transporter superlamfly. J Biol Chem 
272: 10303-10310. 

27. Shroyer NF, Lewis RA, AUikmets R, Singh N, Dean M, et al. (1999) The rod 
photoreceptor .ATP-binding cassette transporter gene, ABCR, and retinal 
disease: from monogenic to multifactorial. Vision Res 39: 2537 — 2544. 

28. Pcelman F, Labeiu C. Vanloo B, Roosbcek S, Dcvaud G, et al. (2003) 
C^haracterization of the .ABCJA transporter subfamily: identitication of 
prokaryotic and eukar^'otic members, phylogeny and topology. J Mol Biol 
325: 259-274. 

29. Zhong M, Molday LL, Molday RS (2009) Role of the C terminus of die 
photoreceptor ABCA4 tiansporter in protein folding, function, and retinal 
degenerative diseases. J Biol Chem 284: 3640—3649. 

30. WengJ, Mata NL, Azarian SM, Tzckov RT, Birch DG, et al. (1999) Insights 
into the function of Rim protein in photoreceptors and etiology of Stargardt's 
disease from the phenotype in abcr knockout mice. CJell 98: 13 — 23. 

31. Quazi F, Lenevich S, Molday RS :2012i .ABC,A4 is an N-retinylidene- 
phosphatidylcthanolaminc and phosphatidylethanolamine importer. Nat Com- 
mun 3: 925. 

32. Sheflield VC, Stone EM (2011) Genomics and die eye. N Engl J Med 364: 
1932-1942. 

33. Mata NL, Weng J, Travis GH (2000) Biosynthesis of a major lipofuscin 
fluorophore in mice and humans with ABCR-mediated retinal and macular 
degeneration. Proc Nad Acad Sci U S A 97: 7154-7159. 

34. Mata NL, Tzekov RT, Liu X, WengJ, Birch DG, et al. (2001) Delayed dark- 
adaptation and lipoluscin accumulation in abcr+/— mice: implications for 



PLOS ONE I www.plosone.org 



7 



March 2014 | Volume 9 | Issue 3 | e91962 



Novel Mutation in ABCA4 in Stargardt Disease 



invoK cmcnt of ABCR in age-related macular degeneration. Invest Ophthalmol 
Vis Sci 42: 1685-1690. 

35. Sparrow JR, Boulton M (2005) RPE lipoluscin and its role in retinal 
pathobiology. Exp Eye Res 80: 595—606. 

36. Wu Y, Fishkin NE, Pande A, Pande J, Sparrow JR (2009) Novel Upofuscin 
bisretinoids prominent in human retina and in a model of recessive Stargardt 
disease. J Biol Chem 284: 20155-20166. 

37. Kim SR Jang YP, Jockusch S, Fishkin NE, Turro NJ, et al. (2007) The aU-trans- 
rctinal dimcr series of lipofuscin pigments in retinal pigment epithelial cells in a 
recessive Stargardt disease model. Proc Natl Acad Sci USA 104: 
19273-19278. 

38. Ben-Shabat S, Parish C^A, VoUmer HR, hagaki Y, Eishkin N, et al. (2002) 
Biosynthctic studies of A2E, a major fluorophore of retinal pigment epithelial 
lipofuscin. J Biol Chem 277: 7183-7190. 

39. Zemant J, Schubert C, Im KM, Burke T, Brown CM, et al. (201 1) Analysis of 
the ABCA4 gene by next-generation sequencing. Invest Ophthalmol Vis Sci 52: 
8479-8487. 

40. Michaelides M, Chen LL, Brandey MAJr, Andorf JL, Isaak EM, et al. (2007) 
ABCA4 mutations and discordant ABCA4 alleles in patients and siblings with 
bull's-eye maeulopathy. BrJ Ophthalmol 91: 1650—1655. 

41. Klevering BJ, Deutman AE, Maugeri A, Crcmcrs FP, Hoyng CB (2005) The 
spectrum of retinal phenotypes caused by mutations in the .\BCA4 gene. 
Graefes Arch CJlin Exp Ophthalmol 243: 90-100. 

42. Cremers EP, van de Pol DJ, \'an Driel M, den Hollander Al, van Haren IJ, et al. 
(1998) Autosomal reeessix e retinitis pigmentosa and eone-rod dystrophy caused 
by splice site mutations in the Stargardt's disease gene ABCR. Hum Mol Genet 
7: 355-362. 



43. Birch DG, Peters AY, Locke KL, Spencer R, Mcgarity CF, et al. (2001) Visual 
function in patients with cone-rod dystrophy (CRD) associated with mutations in 
die ABCA4(ABCR) gene. Exp Eye Res 73: 877-886. 

44. Fishman GA, Stone EM, Eliason DA, Taylor CM, Lindeman M, et al. (2003) 
ABCA4 gene sequence variations in patients with autosomal recessive cone-rod 
dystrophy. Arch Ophthahnol 121: 851-855. 

45. Klevering BJ, van Dricl M, van de Pol D]. Pinrkers AJ, Crcmcrs FP, et al. (1999) 
Phenolypir variations in a family wilh relinal dvstrophy as result of different 
mutations in the ABCR gene. BrJ Ophthalmol 83: 914—918. 

46. Martinez-Mir A, Paloma E, ,\llikmets R, Ajuso C, del Rio T, et al. (1998) 
Retinitis pigmentosa caused by a homozygous mutation in the Stargardt disease 
gene ABCR. Nat Genet 18: 11-12. 

47. Allikmets R, Shroyer NF, Singh N, Seddon JM, Lewis RA, et al. (1997) 
Mutation of the Stargardt disease gene (ABCR) in age-related macular 
degeneration. Science 277: 1805—1807. 

48. Strom SP, Gao YQ, Martinez A, Ortube C, Chen Z, et al. (2012) Molecular 
diagnosis of putative Stargardt Disease probands by exome sequencing. BMC 
Med Genet 13: 67. 

49. Webster AR, Heon E, Lotcry AJ, Vandenburgh K, Casavant TL, et al. (2001) 
An analysis of allelic variation in the ABCA4 gene. Invest Ophthalmol Vis Sci 
42: 1179-1189. 

50. Walker JE, Saraste M, Runswick MJ, Gay NJ (1982) Distandy related sequences 
in the alpha- and beta-subunits of ATP synthase, myosin, kinases and other 
ATP-requiring enzymes and a common nucleotide binding fold. EMBO J 1 ; 

945-951. 

51. Walker JE, Saraste M, Runswick MJ, Gay ^.J (1982) Distantly related sequences 
in the alpha- and beta-subunits of ATP synthase, myosin, kinases and other 
ATP-requiring enzymes and a common nucleotide binding fold. EMBO J 1 ; 
945-951. 



PLCS ONE I www.plosone.org 



8 



March 2014 | Volume 9 | Issue 3 | e91962 



